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ABSTRACT 

In a previous paper we have described a technique to derive constraints on 
the differential emission measure (DEM) distribution, a measure of the temper¬ 
ature distribution, of collisionally ionized hot plasmas from their X-ray emission 
line spectra. We apply this technique to the ChandrajYIETG spectra of all of 
the nine hot stars available to us at the time this project was initiated. We find 
that DEM distributions of six of the seven O stars in our sample are very similar 
but that 6^ Ori C has an X-ray spectrum characterized by higher temperatures. 
The DEM distributions of both of B stars in our sample have lower magnitudes 
than those of the O stars and one, r Sco, is characterized by higher tempera¬ 
tures than the other, /3 Cm. These results confirm previous work in which high 
temperatures have been found for 6^ Ori C and r Sco and taken as evidence for 
channeling of the wind in magnetic fields, the existence of which are related to 
the stars’ youth. Our results demonstrate the utility of our method for deriving 
temperature information for large samples of X-ray emission line spectra. 


1. Introduction 

In hot stars. X-rays are emitted by gas heated in shocks in the wind. In most cases, these 
shocks are thought to result from an instability inherent in the mechanism which drives the 
wind: absorption of stellar UV radiation in line transitions. In some cases, however, other 
physical processes may play a role in the heating of gas to X-ray emitting temperatures. 
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For example, in very young hot stars, channeling of the wind by magnetic helds may play a 
strong role in heating the wind (e.g., ud-Doula & Owocki 2002). Among the evidence that 
X-ray emission from some hot stars may be due to a mechanism other than the instability 
of the radiative driving mechanism, is the fact that the X-ray spectra of some hot stars are 
characteristic of emission by plasma of higher temperature than others. 

Line emission from hot plasmas is highly sensitive to temperature. Therefore, insofar as 
theories of hot star winds make predictions for the temperature of X-ray emitting plasmas in 
them, high-resolution X-ray spectroscopy offers a powerful method with which to test those 
theories. However, there are two obstacles in using high-resolution spectroscopy to study 
the winds of hot stars. First, in most hot stars, the winds from which the X-ray emission 
arises have velocities in excess of 1000kms“^. Owing to this fact. X-ray emission lines are 
broadened and blended, making it difficult or impossible to measure fluxes of many emission 
lines (e.g., Cassinelli et al. 2001; Kahn et ah 2001; Schulz et al. 2003). Second, the theory of 
hot star winds is not sufficiently developed to the point where detailed predictions for X-ray 
emission line spectra can be made for direct comparison with observations. 

While there are open questions regarding the origin of the X-ray emitting plasma in hot 
stars, it is well established that the emitting plasma in hot stars is in collisional ionization 
equilibrium (CIE, see, e.g., Paerels & Kahn 2003). In a previous paper (Wojdowski & 
Schulz 2004, hereafter Paper I), we have described a method for measuring the differential 
emission measure (DEM) distribution, a measure of the temperature distribution of a hot, 
optically thin plasma in collisional ionization equilibrium (CIE). In this method, a number 
of templates, each of which contains all of the emission lines of a single ion, are ht to an 
observed high-resolution spectrum. The best-fit normalization for each of these templates 
then provides a measure of the product of the elemental abundance and a weighted average 
of the DEM distribution over the range where the ion emits. These measurements may 
then be plotted in a way that may be described as a one-dimensional “image” of the DEM 
distribution. This method is well suited to the study of hot stars because the use of templates 
makes use of the information available from line blends and because it does not require that 
assumptions about the form of the DEM distribution be made. 

In §2 we describe our sample of nine early-type stars which have been observed with 
C/iandra/HETGS, the data, and the application of the method of Paper I to it. In §3 we 
discuss the implications of our results. 
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2. Spectral Data Analysis 

We obtained the data from the observations of all of the nine bright O stars observed 
with the C'/iandra High-Energy Transmission Grating Spectrometer (HETGS). In Table 1, we 
list these nine stars, the distances and interstellar absorption columns to them that we adopt, 
and the effective exposure times, Chandra Observation ID numbers, and any references for 
the Chandra/HETGS observations of them. 

For all of the stars in our sample, our analysis proceeds from binned count spectra 
(PHA hies), instrument effective areas and exposure times (ARE hies), and wavelength 
redistribution functions (RMF hies). For 6^ Ori G we use the spectral data and instrument 
response hies described by Schulz et ah (2003). For all of the other stars, we use the 
following procedure. We use spectra extracted as in the standard pipeline processing done 
by the Chandra X-ray Genter. We calculate the appropriate ehective area and exposure 
time for the positive and negative hrst orders of the the High Energy Grating (HEG) and 
the Medium Energy Grating (MEG) for each observation using a standard procedure^. For 
each of the two gratings we add the spectra and ehective areas for the positive and negative 
hrst orders. For those stars observed in two pointings, we add the spectra and exposure 
times and average the ehective areas for the two pointings. In all cases we use the standard 
redistribution functions. The background expected for these spectra is very low and we make 
no correction for it. 

For each of the nine data sets, we apply the procedure described in Paper I — i.e., we 
ht the spectra using templates, each of which contains all of the emission lines of an ion. 
In each of our hts to the stellar data, we use two bremsstrahlung components except in the 
case of r Sco, 9^ Ori G, and ( Ori where three are necessary to ht the continuum. In each 
case, the model spectrum includes neutral, interstellar absorption hxed at the value listed in 
Table 1. As explained in Paper I, the normalization parameter of each emission line template 
obtained in the ht implies a constraint on the product of the elemental abundance and a 
weighted average of the DEM distribution: 

p _ , f ^^»zATmT)d\ogT 
r„t>zJT)dlogT 

where Z is the atomic number of the element, is the charge state, A is the elemental 
abundance dehned relative to solar, V is the diherential emission measure distribution, T is 
temperature, and 6 is the sum of the power functions (described in Paper I) for the lines 
of the ion. In Figures 1 and 2 we plot D vs. Tp, where Tp^z,z is the temperature where 
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Table 1. Stars and Data in our Study 


HR# 

HD # 

Name 

Type 

Adopted Values 

Exposure 

(sec) 

Obs IDs 

References'^ 



Distance 

(pc) 

Column Density 
(lO^'^ cm“^) 

Type Distance Absorb 

Data^ 

3165 

66811 

C Pup 

04I(n)ff 

450 

1.0 

68,598 

640 

24(22) 

10,11 

15 

1 

1895 

37022 

e* Ori O 

04-6p var 

450 

19 

82,975 

3,4 

26^2) 

15 

5 

5,6 

8281*' 

206267A 


06.5V((f))'* 

800 

30 

72,557 

1888,1889 

25(22) 

12 

16 

27 

2782 

57061 

T CMa 

0911 

1480 

or 

bo 

n 

87,095 

2525,2526 

24^2) 

7 

7 

27 

1899 

37043 

t Ori 

09IIP 

440 

2.0 

49,917 

599,2420 

25^2) 

13 

15 


1852 

36486 

6 Ori A 

09.5ie 

501 

1.5 

49,045 

639 

24(22) 

15 

15 

4 

1948 

37742 

C Ori A 

09.7Ibs 

501 

3.0 

59,640 

610 

24^2) 

15 

15 

2 

6165 

149438 

T ScO 

B0.2V 

132 

2.7 

59,630 

638 

23 

8(3) 

15 

3 

4853 

111123 

p Cru 

B0.5IIP 

110 

1.7 

74,379 

2575 

18 

8(9) 

14 



References. — (1) Cassinelli et al. 2001; (2) Waldron Sz, Cassinelli 2001; (3) Cohen et al. 2003; (4) Miller et al. 2002; (5) Schulz 
et al. 2000; (6) Schulz et al. 2001, 2003; (7) Moitinho et al. 2001; (8) Perryman et al. 1997; (9) Alcala et al. 2002; (10) Brandt et al. 
1971; (11) Schaerer et al. 1997; (12) Simonson 1968; (13) consistent with Hillenbrand 1997 and references therein; (14) Code et al. 
1976; (15) Berghofer et al. 1996, BSC; (16) Shull & van Steenberg 1985; (18) Hiltner et al. 1969; (22) Mafz-Apellaniz et al. 2004; 
(23) Walborn 1971; (24) Walborn 1972; (25) Walborn 1973; (26) Walborn 1981; (27) Schulz 2003 

^Previous publications describing these data from Chandra. 

^The designation HR 8281/HD 206267 includes four stars which we resolve. The spectrum is for HD 206267A. 

^derived using A^h = 5.8 x 10^^Eb-v cm~^ 

•^HD 206267A consists of stars of type 06.5V((f)) and BOV in a 3.7 day orbit and a third component of type 08V (Harvin et 
ah, in prep.). 

^spectroscopic binary 

^variable of BY Dra type 

Semission-line star 

^eclipsing binary 

Variable star of beta Cep type 

V determination of type 07V for 6^ Ori C has also been made by Conti (1972) 

'^References in parenthesis indicate a secondary reference. 


4 :^ 
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the function Oz^z peaks. In all of these plots, we use the same temperature range on the 
horizontal axis as in the plots in Figures 2 and 3 of Paper I and, also as in those hgures, the 
vertical axes span three orders of magnitude. Therefore slopes of lines in these various plots 
may be compared. As in Paper I, we model the emission lines as Gaussians and allow the 
centers and widths of all of the Gaussians to vary together. The parameters which describe 
the centers and widths of the lines are, respectively, and Vt and we tabnlate the best-£t 
valnes of these parameters in Table 2. The exact dependence of the model line profiles on 
these parameters is given in equations 16 and 17 of Paper I. While we have not qnantitatively 
assessed the qnality of the fits, it may be seen from the spectral plots that the fits are qnite 
good considering the large nnmber of lines apparent in the data. 

As we have described in Paper I, some excited states of some ions are metastable and 
may undergo further excitation before decaying to ground (pnmping) and this leads the 
Inminosities of some emission lines to depend on physical variables other than the DEM 
distribntion. In the vicinity of a hot star, ions in many metastable states are susceptible 
to excitation by the absorption of UV radiation from the stellar photosphere. As we have 
also described in Paper I, we acconnt for any line pnmping by letting the line ratios in the 
templates vary as with collisional pnmping with one free density parameter for each ion. This 
approach is somewhat ad hoc. However, the fact that we obtain good fits to the spectra 
indicates that pumping is well acconnted for. Therefore, we do not expect significant errors 
in onr constraints dne to line pnmping or inaccnracies in our treatment of it. 

In the plots of D vs. Tp, for nearly every two ions with close values of Tp, the corre¬ 
sponding valnes of D generally appear to be consistent with each other or differ by no more 
than a factor of two, indicating that the emission line spectrnm is consistent with emission 
from a plasma with solar abnndances. One exception to this is the fact that, for ( Pup, the 
values of D for oxygen are approximately an order of magnitnde less than the valnes for ions 
with similar valnes of Tp. A similar resnlt has been obtained by Kahn et al. (2001). The 
valnes of D for nitrogen appear to fit better with the rest of the data points than do the 
valnes for oxygen. Therefore, it is probably the case that oxygen is underabnndant in the 
wind of C Pnp. However, from this analysis alone, it is impossible to exclnde a hypothesis in 
which the DEM distribntion of ( Pup peaks at approximately 3 x 10® K, the abundance of 
nitrogen is greater than the solar valne, and the abnndance of oxygen is near the solar value. 
In addition, there is significant scatter in the data points for C Pnp in the temperatnre range 
4 X 10®-10^ K. Some, bnt not all, of this scatter might be acconnted for by an nnderabnndace 
of iron. For r Sco, like C Pup, the values of D for oxygen fall below the values for other ions 
at nearby temperatnres, thongh not by as much. 

Another exception to onr results being consistent with solar abnndance is 9^ Ori G. For 
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Fig. 1.— Plots of D as a function of temperature for ( Pup, 6^ Ori C, HD 206267A, r CMa, 
L Ori, and r Sco. The horizontal location of the data points is at the temperature where the 
ion’s emission peaks. The horizontal extents of the diamonds give the temperature range 
over which the ions’ emission is significant and the vertical extents of the diamonds give the 
statistical conhdence intervals for the values of D. The designation “(ul)” indicates an upper 
limit. 
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Fig. 2.— DEM constraints for C Ori A, 6 Ori, and /3 Crn. 
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Fig. 3.— MEG spectral data (black) and best-fit models (red) for ( Pup, 9^ Ori C, 
HD 206267A, r CMa, l Ori, and r Sco. We include a number of line labels though not 
all of the labeled lines are detected in every spectrum. 
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Fig. 4.— MEG spectral data and best-fit models for ( Ori A, 5 Ori, and Cm. 
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this star, most of the values of D appear to lie on a single curve except for those of the 
ions Fe XX- XXIV which lie below it by a factor of approximately four. Schulz et ah (2003) 
£t this spectrum using an explicit parameterized empirical model for the DEM distribution 
with variable element abundances. The best-fit DEM distribution these authors derive has 
three peaks at 7.9, 25, and 63 MK and troughs at 13 and 40 MK and vanishes towards 
temperatures higher and lower than the three peaks. This three-peaked structure is not 
immediately apparent in our constraints for 9^ Ori C shown in Figure 1. We have explicitly 
evaluated the right-hand side of equation 1 for the element abundances and DEM distribution 
of Schulz et ah (2003)^ and while the values of D for some ions have statistically significant 
differences, there does not appear to be an overall systematic difference between the values 
of D obtained by the two methods. In particular, as with the values of D obtained by our 
method, the triply peaked DEM distribution is also not immediately apparent in the values 
of D obtained from the abundances and triply peaked DEM distribution of Schulz et al. 
(2003). This is somewhat puzzling. However, while a triply peaked DEM distribution is not 
immediately apparent in our results, from a closer look it may be seen how a triply peaked 
DEM distribution may be found with the method of Schulz et al. (2003). In the method of 
Schulz et al. (2003), the formula for the quality of fit includes a penalty for variations in the 
DEM distribution so, if possible, elemental abundances will be adjusted to produce a flat 
DEM distribution. Variations in the DEM distribution will be found only when lines of the 
same element require it. Eollowing our values of D for Ee XVII-Fe XXIV, the first peak, the 
following trough, and the rise to the second peak may be seen. The trough after the second 
peak may be seen in the drop from Ar XVII to Ar XVIII and from Ca XIX to Ca XX. The 
rise to the third peak may be seen from the fact that the point for Ee XXV is higher than 
the one for Fe XXIV. 

In order to compare our constraints on the DEM distributions of the stars in our sample 
with each other, we plot a composite of the nine stars in Figure 5. This figure is constructed 
as follows. For each star, we ignore those values of D without lower limits. For the case of 
C Pup we also ignore the values of D for oxygen. Then, for those ions with the same value 
of Tp as given in Table I of Paper I, we average the values of D weighting by the inverse 
square of the uncertainties. The lines connect the averaged values of D for individual stars. 
Because the method used to construct this figure is ad hoc and because doing so would make 
the hgure harder to read, we do not attempt to include errors in this figure. However, the 
magnitude of the errors for this hgure may be inferred from Figures I and 2 and may also be 


^Schulz et al. (2003), in their Figure 7 plot the “DEM”. However, this is not the DEM as we define it 
(emission measure per temperature decade) but the emission measure per temperature bin (with A log T = 
0.05). Therefore, it must be multiplied a factor of 20 to be compared with our results. 
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inferred from the jagged curves. The lines for HD 206267A, l Ori, C, Pup, C, Ori A, r CMa, 
and 5 Ori A are very similar and difficult to show together. Therefore, we show the line for 
C Pup and a hlled polygon that encloses the lines for the other hve stars. The axes and aspect 
ratio of this plot are chosen so that slopes in this hgure can be compared with Figures 2 and 
3 of Paper I and Figures 1 and 2 of this paper. With the exceptions of r Sco and 9^ Ori C, 
it may be inferred from these plots that, in the temperature range 10®-10®'^K, the DEM 
distributions have peaks at temperatures of a few times 10® K and decrease towards higher 
temperatures. The DEM distribution of (3 Cm, in addition to having a smaller magnitude 
than the other stars, declines more rapidly towards high temperatures. For r Sco, the peak 
of the DEM distribution and the majority of the emission appear to be at approximately 
10^ K. Most of the emission measure for 9^ Ori C appears to he near 3 x 10^ K. In Figures 3 
and 4, it may be seen that for many of stars, the points for some or all of the ions Fe XXIV, 
Fe XXV, and Fe XXVI are higher than those of lower charge states of iron. This suggests 
the possibility that these DEM distributions increase at temperatures near 10® K. However, 
all of these high temperature points have low statistical signihcance and it is not possible 
from this analysis to conhrm the presence of any such high temperature increase in the DEM 
distribution. 

We give the best-£t line shifts and widths from our fits in Table 2. If we add the 
systematic uncertainty of the Chandra HETG absolute wavelength calibration (a few tens 
of kms“^)® to the tabulated statistical confidence intervals, all of the line shifts we measure 
are consistent with zero, with the exception of C Pup and the possible exception of C Ori. 
The blueshift we find for ( Pup is consistent with those reported by Kahn et al. (2001) and 
Cassinelli et al. (2001). The latter authors attribute this blueshift to obscuration of the 
far side of wind by photoelectric absorption in nearer parts of the wind. The blueshift we 
observe for C Ori A may also be due to this effect. This possibility is discussed by Waldron 
& Cassinelli (2001). 


3. Discussion 

We use a method, described in Paper I, to determine the DEM distribution for a sample 
of nine hot stars which works by fitting entire X-ray spectra. The fact that for all of the 
nine stars in our sample good spectral fits are obtained suggests that the method is free of 
significant systematic errors. We find that the magnitudes of the DEM distributions of the 
two B stars r Sco and (3 Cm are much lower than those of the O stars as expected from 


^http: / / cxc.harvard.edu/proposer/POG/ 
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Fig. 5.— A composite approximate representation of of the DEM constraints we derive for 
the nine stars. The hlled polygon encloses the lines for HD 206267A, l Ori, Q Pup, C, Ori A, 
r CMa, and 5 Ori A which are very similar and difficult to show together. The jagged shape 
of the curves gives some idea of the magnitude of the error for the DEM constraints. Details 
of the construction of this hgure are given in the text. 


Star 

Type 

Pr (kms 

Vt (kms 

C Pup 

OSIaf 

-460±20 

770±20 

Ori c 

06pe 

_29+^® 

326±18 

HD 206267A 

06.5V((f)) 

-10±140 

1020+1^° 

C Ori A 

09 lab 

-130±20 

689±19 

r CMa 

09Ib 


990^^°° 

i Ori 

09HI 

60±40 

930^^° 

5 Ori A 

09.511 

10±30 

690±30 

r Sco 

BOV 

37±10 

159^}^ 

f3 Cm 

B0.5HI 

-12±19 

190±30 


Table 2: Fit Line Shifts and Width Parameters 


Note. — v-c and are, respectively, the shift and Gaussian a, see Paper I. 
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previous results (e.g., Berghofer et al. 1996). For all of the stars in our sample except for 
6^ Ori C and r Sco we find the DEM distributions to be peaked at a temperature of a few 
times 10® K or less, to decline up to several times 10^ K. Though we cannot conhrm it, the 
DEM distributions of these stars may also increase at temperatures near 10® K. In contrast, 
we hnd that the peak of the DEM distribution of r Sco is near 10^ K and the peak of the 
DEM distribution of 6^ Ori C is near 3 x lO'^ K. For those two stars we also hnd that the 
emission lines are narrow compared to stars of similar stellar type. High temperatures and 
narrow lines have previously been found for r Sco by Cohen et al. (2003) and for 6^ Ori C 
by Schulz et al. (2001, 2003) using these same data from Chandra. In both cases, it has 
been noted that the stars are young and that the high temperatures and narrow lines have 
been attributed to channeling of the wind in fossil magnetic helds. In this scenario, wind 
travels up the sides of magnetic loops, then shocks as wind from the two sides of the loop 
collides. Because the post-shock gas is stationary, narrow lines are predicted. The fact that 
the post-shock gas is stationary also resnlts in a larger velocity difference between the gas 
before and after the shock and, therefore, higher temperatures (ud-Doula & Owocki 2002). 

The fact that the X-ray emitting plasma in the wind of r Sco is hotter than the X- 
ray emitting plasma in most stars has been apparent from spectra of much lower resolution 
(e.g., Cohen, Cassinelli, & MacFarlane 1997a; Cohen, Cassinelli, & Waldron 1997b). A 
hot X-ray emitting plasma has also previously been detected from the Orion Trapezium 
cluster which contains 6^ Ori C (Yamauchi et al. 1996) thongh it was not nntil a high- 
resolntion image was obtained with Chandra by Schulz et al. (2001) that this hot plasma 
could be conclusively associated with 9^ Ori C. Therefore, our Ending that the temperatures 
of the X-ray emitting plasmas 9^ Ori C and r Sco are high is neither a new result of high- 
resolution spectroscopy nor of our analysis techniqne. However, as we have shown here, the 
method described in Paper I is useful in deriving well-defined qnantitative constraints on 
elemental abundances and differential emission measnre distributions from high-resolution 
X-ray spectra and, especially, large samples of high-resolution X-ray spectra. We expect 
that these results and the results of similar observations and analyses will provide direction 
for developing theory on X-ray emission in hot stars. 


We thank John Houck for assistance implementing our analysis technique in ISIS, Dan 
Dewey for a careful reading of the manuscript, David Huenemoerder for providing data for 
the comparison of our results for 9^ Ori C with those of Schulz et al. (2003) and the referee, 
Ehud Behar, for helpfnl comments. This research has made use of the National Aeronautics 
and Space Administration’s (NASA’s) Astrophysics Data System Bibliographic Services, the 
SIMBAD database, operated at CDS, Strasbonrg, France, and the Chandra Data Archive, 
part of the Chandra X-Ray Observatory Science Center (CXC) which is operated for and on 
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